Edge mass current and roles of Majorana Fermions in superfluid '^He A-phase 
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The total angular momentum associated with the edge mass current flowing at the boundary in 
the superfluid ''He A-phase confined in a disk is proved to be L = Nh/2, consisting of L^'' = Nh 
from the Majorana quasi-particles (QPs) and L"^""* = —Nh/2 from the continuum state. We show it 
based on an analytic solution of the chiral order parameter for quasi-classical Eilenberger equation. 
Important analytic expressions are obtained for mass current, angular momentum, and density of 
states (DOS). Notably the DOS of the Majorana QPs is exactly No/2 (TVq: normal state DOS) 
responsible for the factor 2 difference between L^'' and L"^""*. The current decreases as E~'^ against 
the energy E, and L{T) oc — T^. This analytic solution is fully backed up by numerically solving 
the Eilenberger equation. We touch on the so-called intrinsic angular momentum problem. 

PACS numbers: 67.30.hp, 67.30.ht 
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The superfluid is a paradigmatic and fertile testing 
ground to investigate the central questions concerning 
the origin of anisotropic superfluidity since the p-wave 
pairing symmetry for the ^Hc A-phase and B-phase is 
firmly established [H. Recently the topological aspect of 
the '^He A- and B-phases has become a focus because 
there has been no other concrete topological superfluids 
or superconductors established so far. The topological 
superfluid defined by a topological number 0. At an 
edge of the superfluid ^Hc faced to a topologically triv- 
ial vacuum, a superfluid gap is closed by a topological 
phase transition at the interface. The surface Andreev 
bound state of Majorana Fermions must appear Q . The 
topological features arc quite different in the '^He A- and 
B-phases. The "^He A-phase is a chiral superfluid with the 
spontaneous edge mass current while B-phase is a helical 
superfluid with the edge spin current The Majorana 
nature for the A- and B-phases is also distinctive; Ma- 
jorana Fermions have the surface Andreev bound state 
with linear dispersion relation centered at the zero en- 
ergy forming a "Majorana valley" in the A-phase 0, [Hj 
and a "Majorana cone" in the B-phase This Ma- 

jorana cone is emerging by recent experiments 0, [lo| . 
However, there has been no firm evidence for detecting 
the Majorana zero mode. We are still in the process to 
better characterize this elusive quasi-particlc. 

The edge mass current spontaneously generated at the 
boundary in the superfluid '^He A-phase is carried both 
by the Majorana quasi-particles (QPs) within the energy 
gap and by QPs in the continuum state outside the gap. 
Those two kinds of QPs play a different role for the edge 
mass current and associated total angular momentum L 
of the whole system, such as a disk. Stone and Roy [11 1 
proved that in the A-phase the magnitude of the total 
angular momentum by the edge mass current is Nh/2 
at the zero temperature {N is the total number of '^He 
atoms in disk geometry). 

Here we pose problems on what amount of this mass 
current or the angular momentum L carried by the Ma- 



jorana QPs and the remaining continuum state. It will 
turn out that exactly i^'' — Nh for Majorana QPs and 
^cont _ —Nh/2 for the continuum state, giving rise to 
L — Nh/2 for the total angular momentum of the system. 
Then it is shown that the excess angular momentum due 
to the Majorana QPs comes from the fact that the den- 
sity of states (DOS) of the Majorana QPs equals exactly 
to Nq/2 with the normal state DOS iVg at the Fermi 
energy Ep. It will be also demonstrated that the con- 
tribution of the particles with the energy E to the mass 
current decreases as E~^ , namely it is neither confined 
to the narrow energy shell region around the Fermi level, 
nor the whole particles deep inside the Fermi level. It is 
marginal. 

We base our arguments with the quasi-classical Eilen- 
berger theory that is vahd for ^ > fc^^ (^ - 10-100 
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nm the coherent length and k~p' 
isfied for the superfluid '^He. We first find an analytic 
solution for the order parameter for the A-phase when 
the system has a boundary. This analytic solution satis- 
fies the Eilenberger equation and gives useful and trans- 
parent information on various physical properties at the 
boundary. Since this analytic solution exactly satisfies 
the Eilenberger equation, but is not a self-consistent so- 
lution for the whole problem unless taking the infinite 
cutoff energy, we back it up by numerically solving this 
problem self-consistently, finding that the analytic solu- 
tion captures the essential points and the numerical solu- 
tion gives additional insights to further characterize the 
Majorana QPs. 

We will touch on the so-called intrinsic angular mo- 
mentum (lAM) problem, a long standing controversy 
concerning L at the zero temperature. They are sum- 
marized as L^*^ = {Nh/2M/SJEFV with 7 = 6 [iMi, 
7 = 1 and 7 = 2 [llliB], where A is a superfluid 

gap. Since /\/Ef ~ lO'^ for the superfluid ^He, the lAM 
can be observed macroscopicaUy only when 7 = 0. It is 
our hope that the present investigations may contribute 
to this lAM problem. 
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We consider disk geometry with a radius R much larger 
than the coherence length and a small thickness along 
the z-direction to align ^-vector toward the z-direction, 
where c?-vector is known to be also aligned toward the 
z-direction by the dipole interaction This system 
is realized by confining the superfluid '^He A-phase in 
a slab of sub-/xm thickness Q that is already realized 
experimentally (2ll |. Since the d- vector is fixed, we con- 
sider the orbital part of the order parameter A(a;, k) = 
Ax{x)kx + iAy{x)ky in the one-dimension toward the ra- 
dial direction x, where k normalized by kp is the relative 
momentum of a Cooper pair. The coefficients and 
Ay can be chosen as a real number without loss of gen- 
erality. We assume that the edge at a; = is specular 
where only the coefficient A^ is suppressed at the edge 
and \ Ax\ = \Ay \ far from the edge because the chiral state 
is recovered in the bulk. 

Microscopic information of the edge mass current 
is contained by the quasi-classical Green's functions 
g{x,k,ujn), /(x, fc, a;,i), and f{x,k,ujn) which satisfy the 
normalization condition = 1 — //. The quasi-classical 
Green's functions are calculated by using Eilenberger 
equation [22| as 



ujn + hvpkx-^^ f = A{x,k)g, 



(1) 



where vp is the Fermi velocity and a;„ = (2n+ VjirkBT is 
the Matsubara frequency with n G Z. The self-consistent 
condition for the pair potential is given as the gap equa- 
tion 

A(a;,fc) = NonkeT 

X {V{k,k')[f{x,k\uJn)+l{x,k',OJnr])^,, 



0<aj„<aJc 



(2) 



where {• ■ •)k indicates the Fermi surface average and Wc is 
a cutoff energy set to be ujc = 407rfcBTc in later numerics 
with the transition temperature T^. The pairing inter- 
action is given as V{k, k') = 3gik ■ k' for Cooper pairs 
with the orbital angular momentum / = 1, where gi is a 
coupling constant for the p-wave channel. We use the re- 
lation {giNa)-^ = MT/Tc) + 2TrkBTy]a^„<^,^n^ and 
solve Eq. ^ by the Riccati method |23|.124| in numerical 
calculations. 

By using the quasi-classical Green's function, the mass 
current is calculated by 



j{x) = mvpNoirkBT 



E 



{kg{x,k,ujn))k, (3) 



where m is the mass of the ■^He atom. The mass current 
and local density of states (LDOS) for an energy E are 



given by 

jix,E) = {jix,k,E))k 

= mvpNoikRc [g(x,k,ujn)\i^^^^E+ir,])k, (4) 
N{x,E) = {N{x,k,E))k 

= NQ{Re[g{x,k,ujn)\i^^^^E+ir,])k, (5) 

respectively, where 77 is a positive infinitesimally small 
constant. 

In order to better characterize the Majorana QPs, it 
is essential to decompose the edge mass current into the 
Majorana QPs of the Andrcev bound state within the en- 
ergy gap and the continuum state outside the gap. This 
is our first task. 

We solve Eilenberger equation ([!]) under the pair po- 
tential with 



A^ = Ao tanh [ - ] , Ay = An 



(6) 



where Aq is the superfluid energy gap in the bulk and 
the coherent length is deflned by ^ = hvp /Aq. This 
pair potential is hinted by Ref. 25. This order parameter 
form embodies the fact that the k^- (fcy-)component is 
suppressed (intact) at the boundary. The quasi-classical 
Green's function 1261 



g{x,k,uJn) 



1 



+ Al sin^ 



Aq sin^ 6 cos^ cj) 
2(w„ -|- iAq sin 9 sin ( 



sech'^ ( ^ 



(7) 



is shown to satisfy Eq. ([T|) after lengthy but simple cal- 
culations, where k^ = sin 9 cos (f) and ky = sin sine/) in 
the spherical coordinates on the unit Fermi surface. By 
the quasi-classical Green's function in Ec^. we can 
calculate 6'-anglc-resolved LDOS by Eq. ^ as 



N{x,9,E) 



2-K 



N{x,k,E) 



— - scch^ 
2 



(8) 



for the Majorana bound state \E\ < Aq sin 6* and 
N{x,9,E)^No 




£2 _ A2 sin^ I 



\E\ 



£;2 _ A2 sin^ 9 



scch^ ( ^ 



(9) 



for the continuum state \E\ > Aosin6'. Thus, QPs feel 
the pair potential Ao sin 9 where 9 is the polar angle. The 
6'-angle-resolved LDOS at the edge 2; = for 61 = 7r/2 is 
shown in Fig. [ija). The LDOS from the Majorana zero 
energy mode has a constant and finite value Nq/2. Simi- 
lar value is also obtained by self-consistent calculation in 
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a finite temperature [3, HI- This means that the linear dis- 
persion of the Majorana QPs is constant and independent 
of 9. According to the numerical calculation , the ve- 
locity wmj of the Majorana QPs is vm^ ~ Aq /hkp <C wf- 
We can also calculate 0-angle-resolved mass current 
along the edge by Eq. as 



jy{x,0,E) 



2tt 

mvpNo E 



jy{x, k,E) 



Ao 



scch 



(10) 



for the bound state \E\ < Aq s'mO and 



jy{x,9,E) 



mvpNo 1 E 
4 A^Wl 
E^ 



E^ - A2 sin2 I 



E'^ - Al sin^ 



2\E\ 



sech 



(11) 

for the continuum state \E\ > Aosinfl. 

The 0-angle-rcsolvcd mass current at the edge for 9 ~ 
7r/2 is shown in Fig. [TJb). The Majorana contribution 
to the current is linear in E. The edge mass current in 
the continuum state decreases away from |Ao sin0|. The 
asymptotic behavior of the edge current (11) at a; = is 
estimated as 



jy{x^0,9,E) 



mvpNQ 
16~ 



Aq sin ( 
E~ 



(12) 



for E — Aq sin 6*, implying that the contribution of the 
edge current decreases ~ This power law behavior 

of the decrease means that it is neither confined near Ep, 
nor deep up to the band bottom, rather it is marginal. 

Since QPs fill the energy state up to the Fermi energy 
at the zero temperature, the mass current along the edge 
from the Majorana bound state is obtained analytically 
as 



dEjy{x, k,E) 

Ao sin 9 / /g 

mvpNoAo . 2 / 
6 U, 



(13) 



and that from the continuum state is 

/ r — Ao sin 9 



jr(^) 



dEjy{x, k, E) 



''''' u 



(14) 



where a similar result to Eq. (jl3p is obtained in connec- 
tion with the chiral superconductor Sr2Ru04 27[. Thus 
those currents flow oppositely. In the disk with R ^ 
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FIG. 1. (Color online) Energy profiles of S-angle-resolved 
LDOS (a) and mass current along the edge (b) at a; = 
for Q = 7r/2. At the zero temperature, quasi-particles fill the 
colored (shaded) states in (a). Mass current from the bound 
and continuum states is derived by integrating the blue (light 
gray) and pink (gray) regions in (b), respectively. 



the angular momentum by the edge mass current from 
each state is calculated as 



-MJ 



Nh 



(15) 



The factor 2 difference in those contributions may be 
related to our finding that the DOS of the Majorana QPs 
is Nq/2. Finally, the total angular momentum becomes 
simply to 



Nh 



(16) 



which corresponds to the intrinsic angular momentum 
for 7 = mentioned above. The total angular momen- 
tum due to the edge mass current coincides with that in 
Ref. 11. Here we emphasize that we decompose it analyt- 
ically into two contributions: the Majorana bound state 
and continuum state. Interestingly, a half of the angular 
momentum from the Majorana bound state is canceled 
partly by that from the continuum state, resulting in 
the expected angular momentum. Note that the total 
number of "^He atoms N emerges from the normal DOS 
A^o ~ (3/TOw|-)(7V/y), where V is the volume of the disk. 

Next, we solve Eilenbergcr equation ([1]) and gap equa- 
tion ^ self-consistently and compare to the solution un- 
der the uniform pair potential. The uniform pair poten- 
tial means A(x,k) = A^ikx + iky), where amplitude of 
Aq is determined by the gap equation. In Fig. [21 the 
profiles of the edge mass current at T = O.STc under self- 
consistent pair potential (solid line) and uniform pair po- 
tential (dashed line) calculated by Eq. ^ arc shown. The 
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FIG. 2. (Color online) Profiles of mass current along the 
edge at T = O.STc under self-consistent pair potential (solid 
line) and uniform pair potential (dashed line). The units are 
^0 = hvp /2nkBTc and jo = mvpNonkBTc. Inset: A profile 
of the self-consistent pair potential. 
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FIG. 3. (Color online) Left axis: Temperature dependence of 
angular momentum under self-consistent pair potential (solid 
circles) and uniform pair potential (open circles) . Right axis: 
Temperature dependence of the component of the superfluid 
density tensor p^Jn (solid line) and Pg± (dotted line). 



profile of the edge mass current under the self-consistent 
pair potential is varied gradually over the region within 
^ 10^0 owing to the variation of the pair potential (the 
inset of Fig. [2]) . In contrast the edge mass current under 
the uniform pair potential is localized near the edge. A 
profile of the self-consistent pair potential is also shown 
in the inset of Fig. [2] Because of the specular boundary 
condition, the fca;-component becomes zero at the edge. 
The fcy-componcnt is enhanced by compensating for the 
loss of the fc^-componcnt on the edge, where the polar 
state is realized. In other words, the inverse chiral state 
—kx + iky is mixed with the chiral state + iky near the 
edge. This inverse chiral state must have the Att phase 
winding to minimize the free energy in the axisymmetric 
disk 

The temperature dependence of the angular momen- 
tum Lz{T) by the edge mass current under the above 
self-consistent and the uniform pair potential is shown 
in Fig. |3] by circles. The component of the superfluid 
density tensor parallel (perpendicular) to the direction 
of the point nodes (p°_l) [l^l is also depicted by two 
lines. In both pair potential cases, Lz tends to Nh/2 to- 
ward the low temperature limit and zero at the transition 
temperature. In the uniform case, Lz{T) (open circles) 
has the same temperature dependence as (T) (solid 
line) [2^, because thermal excitations from the point 
nodes reduce the superfluid density. However, Lz{T) in 
the self-consistent case (solid circles) is larger than that 
in the uniform case (open circles) in intermediate temper- 
atures. The reason for the increment of Lz{T) is that the 
additional inverse chiral state with the 47r phase winding 
yields the extra current. The different T-dependences 
between L^iT) and p^^^\^(T) should be observed experi- 
mentally in wide intermediate T-region. 

According to the above analytic solution (O, we can 
derive the low temperature behavior of the total angular 



momentum as 



Trk^T 



o 



TTkgT 



(17) 



with /3 = 1. Note that /3 = 2/3 in the absence of 
point nodes, namely in the two-dimensional Fermi sur- 
face model. The T^-dependence comes from the Ma- 
jorana QPs because jy[x,6,E) in the Majorana bound 
state has linear energy dependence (Fig. [Hb)). Thus, 
the observation of the T^-behavior could establish the 
existence of the Majorana QPs. 

We point out another possibility to experimentally ob- 
serve the Majorana QPs contribution L^'^ and the con- 
tinuum contribution separately by controlling the 
surface roughness. The surface roughness can be con- 
trolled experimentally by coating the surface with '''He 
atoms [30[ • The angular momenta L^"^ and 1,'=°"* can be 
suppressed by the different rates as a function of the sur- 
face roughness. Here we notice in the case of the super- 
fluid ''He B-phase that the surface Majorana bound state 
remains at a partially specular surface [l^l ■ The con- 
crete calculation to estimate quantitatively this change 
is a future problem. 

Recently a paper by Sauls [sij has appeared, in which 
the author discusses the same topics as our own. How- 
ever, he uses the two-dimensional Fermi surface model 
under the uniform pair potential. We reach the same con- 
clusions regarding the angular momentum at zero tem- 
perature with the specular boundary. 

In summary, we have found an analytic solution ([7]) of 
Eilenberger equation for the superfluid '^He A-phase at 
the boundary. This analytic solution is useful enough to 
elucidate the elusive Majorana QPs in a realistic situation 
accessible experimentally. Namely, the analytic expres- 
sions for the DOS, the mass current and the associated 
angular momentum of the disk, that make clear the roles 
of the Majorana QPs. We have succeeded in separating 
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L = Nh/2 into L^-^ = Nh and L™"* = -Nn/2 and also 
in proving the DOS of the Majorana QPs is exactly the 
half of the normal DOS Nq/2. The self-consistent nu- 
merical computation is consistent with those conclusions 
and provides further information on the Majorana QPs 
that facilitates future experiment for detecting this exotic 
particle. 

Finally, we mention that the relation between the in- 
trinsic angular momentum and the edge mass current. 
The mass current by a Cooper pair is cancelled by over- 
lapping Cooper pairs in the bulk. The total current is 
zero by the particle- hole symmetry [2^. At the edge, 
however, since the cancellation does not occur com- 
pletely, the angular momentum L = Nh/2 emerges via 
the edge mass current. The superfluid ^He A-phase in 
disk geometry has the net angular momentum L = Nh/2. 
In three-dimensional thick slabs, the net angular momen- 
tum reduces by the canting of /-vector at the edge [1] . 
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knowledges the support of the Research Fellowships of 
the Japan Society for the Promotion of Science for Young 
Scientists. 
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